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Faraday rotation induced B-modes can provide a distinctive signature of primordial magnetic fields 
because of their characteristic frequency dependence and because they are only weakly damped 
on small scales, allowing them to dominate B-modes from other sources. By numerically solving 
the full CMB radiative transport equations, we study the B-mode power spectrum induced by 
stochastic magnetic fields that have significant power on scales smaller than the thickness of the last 
scattering surface. Constraints on the magnetic field energy density and inertial scale are derived 
from WMAP 7-year data, and are stronger than the big bang nucleosynthesis (BBN) bound for a 
range of parameters. Observations of the CMB polarization at smaller angular scales are crucial to 
provide tighter constraints or a detection. 



I. INTRODUCTION 

Many early universe scenarios predict the existence of cosmological magnetic fields and several observational tech- 
niques are currently being employed to detect such fields [1, 2]. There are also recent claims for the detection of an 
inter-galactic magnetic field [3]. With further confirmation and refinement, these observations can be turned into a 
probe of the fundamental interactions that are necessary to generate a primordial magnetic field and to the physics of 
the early universe. A primordial magnetic field can also have important implications for the formation of first stars, 
growth of cosmic structure and the present universe. 

Of particular interest to us are magnetogenesis scenarios based on phase transitions such as the electroweak phase 
transition, when Maxwellian electromagnetism first emerged [4-11]. In this scenario particle physics uncertainties are 
minimal and magnetogenesis is tightly related to the creation of matter, or "baryogenesis" - the magnitude of the 
magnetic helicity density is approximately equal to the baryon number density. Since the cosmic number density of 
baryons is known, the scenario enables a prediction for the magnetic helicity density that is largely independent of the 
details of the electroweak model. Furthermore, the left-handed magnetic helicity is a direct outcome of parity violation 
present in the electroweak model and, like baryon number, is a remarkable large-scale manifestation of a microscopic 
symmetry violation. 

Non-vanishing helicity has important consequences for the evolution of a magnetic field. Although the field is 
generated on small scales, magnetic helicity allows for an "inverse cascade" where power is transferred from small 
to large scales, resulting in magnetic coherence on larger scales. In the case of magnetic fields generated at the 
electroweak scale, the final coherence scale can be on the parsec to kiloparsec scales [6, 12]. While the coherence 
scale is large, it is not as large as for fields that might be generated during an inflationary epoch. For purposes of 
calculating observational signatures, magnetic fields generated at a phase transition are stochastic. An important aim 
of the present work is to find distinctive observational signatures of magnetic fields that are stochastic on sub-Mpc 
coherence scales. 

Once magnetic fields have been injected into the cosmological plasma, the subsequent evolution is described by 
magneto-hydrodynamical (MHD) equations in an expanding spacetime. Power on very small length scales is expected 
to be erased by dissipational mechanisms. Power on very large scales is also suppressed because the magnetic field 
is injected on microscopic scales and cannot extend to arbitrarily large scales. Hence the spectral distribution of the 
magnetic field is expected to decay fast on small and large scales and be peaked on some intermediate scale, which is 
presumably at the parsec to kiloparsec scale. The spectral form of the magnetic field has been investigated recently 
in some detail in Refs. [13, 14] and the results are schematically depicted in Fig. 1. 

On the observational front, magnetic fields within galaxies and clusters of galaxies have been studied for many 
decades and their origin - cosmological versus astrophysical - remains unsettled, though a hybrid explanation is also 
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FIG. 1: The Fourier amplitude b{k) - as given by the square root of the power spectrum - of a stochastic magnetic field 
generated during a cosmic phase transition is expected to grow for k < kj where ki is an inertial scale, and obey some other 
power law for kj < k < k^iss, where fcdiss is a dissipative scale. The work of [13] suggests k^ growth at small k, then k° behavior 
until fcdiss, beyond which the amplitude falls off very quickly. 

conceivable. More recently, observations of TeV gamma ray sources have been used to place lower bounds ~ 10^^® G 
[3, 15-18], and perhaps a measurement ~ 10~^^ G [19], on a magnetic field in inter-galactic space, ~ 10 Mpc away 
from the TeV gamma ray source. It seems likely that such an inter-galactic magnetic field, unassociated with cosmic 
structure, is primordial, but an astrophysical origin, say based on the expulsion of magnetic fields from active sources, 
may also be viable. 

A detection of magnetic fields in the cosmic microwave background (CMB), for example due to Faraday rotation 
(FR) of the CMB polarization, would unambiguously point to a cosmological origin because there are no confounding 
magnetized structures at last scattering. However, an observed FR of the CMB could also be due to magnetic fields 
along the line of sight, especially within the Milky Way. Hence, it is necessary to find distinctive signatures of FR that 
occurred at recombination versus that which happened more recently. As we will see, in addition to its characteristic 
frequency dependence, FR induced B-modes are only weakly damped on small angular scales (high £), which means 
they are likely to dominate B-modes from other sources. 

Earlier work on FR of the CMB in cosmic magnetic fields has largely focused on the effect of a uniform magnetic 
field [20, 21] and, when a stochastic magnetic field has been considered, a "thin" LSS was often assumed [22-24]. This 
approximation results in a tremendous technical simplification but it is not clear if it is suited to study the effects 
of a magnetic field with coherence scale smaller than the thickness of the LSS. In [25], the finite thickness effects 
of the LSS were modelled by approximating the visibility function with a Gaussian profile. As we show, the thin 
LSS approximation is sufficient for order of magnitude estimates, but can be wrong by factors of a few and in an £ 
dependent way. 

In the present paper, we study the effect of a primordial magnetic field on the CMB polarization, focussing on the 
effect of a stochastic field with power on small (sub- Mpc) scales. We calculate the B-mode correlator, C^^, induced 
by a primordial magnetic field. In the thin LSS approximation of Sec. IV, we show that Cj^^ is directly related to 
Cf^ multiplied by the correlation function for the FR, C"", which we calculate in Sec. III. In Sec. V we move on 
to our main calculation of with a thick LSS. Here we find that the physics of FR during recombination can 

be described in terms of "window functions" through which the magnetic field spectrum (see Eq. (65)) appears in 
. We have to resort to extensive numerical efforts to compute the window functions. Our results are described in 
Sec. VI. The window functions themselves are independeirt of the magnetic field power spectrum and are shown in 
Fig. 2. When the window functions are convolved with the magnetic field power spectrum, we obtain Cf^ . In Figs. 3, 
4 we show the results for C^^ for a scale-invariant and a blue magnetic power spectrum. We conclude in Sec. VII. 
We have made our window functions publicly available at http://www.sfu.ca/'-^levon/faraday.html along with a 
short Fortran code for calculating Cj^^ . 



A statistically homogeneous and isotropic stochastic magnetic field is described by the two-point correlator in 
Fourier space as 



II. STOCHASTIC MAGNETIC FIELDS 



(6,(k)6j(k')) = (2^)3^(3) (k -I- k')[(*», - kk)S{k)^ie^jikA{k)] , 
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where S{k) and A{k), the symmetric and anti-symmetric magnetic power spectra, are real functions of = |k| 
Throughout the paper, we use Gaussian CGS units. The energy density in modes up to some value of k is given by 
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The results in [13, 26] suggest the exponents n = 5/2 and n' — 3/2. In our analysis, along with these values, we 
will also consider the case of a nearly scale invariant spectrum, with 2n = 2n' — 0.1, motivated by an inflationary 
mechanism of generation of magnetic fields [27, 28]. 

We will only consider the effect of magnetic field on the CMB for modes with I < Zmax = iC, as computations 
at higher I are very expensive. This corresponds to a minimum comoving scale of 1 Mpc, or fcmax ^ 1 Mpc~^. For 
magnetic fields generated at the electroweak phase transition, the coherence scale is estimated at kpc scales or less 
[6, 12]. Hence fc/ may be expected to be 10'^ Mpc^^. This is much larger than fcmax and thus the I < 10^ modes of the 
CMB are likely to be affected only by the magnetic field modes in the inertial range < k < kj , and the form of the 
power spectrum for k > ki plays no direct role. However, even then, the large k part of the power spectrum would 
still enter when we derive a constraint on the magnetic field, since the amplitude, S**, is fixed by the total energy 
density in the magnetic field. We now make this point clearer. 

First define an "effective magnetic field", -Bcff, in terms of the total energy density in the magnetic field, eo; 



Stt 



(4) 



In other words, B^s is the field strength of a uniform magnetic field that has the same total energy density as our 
stochastic magnetic field. Our constraints will be written in terms of -Bcff- To connect to the amplitude of the power 
spectrum, we first evaluate the energy density in the magnetic field. From Eq. (2) we get 
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The total energy density in the magnetic field, eo, is found by setting k — >■ oo, which is the same as tM{k) for k > kdh 
Thus, we can write 



where 



K=l 
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and hence 
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For a fixed exponent n, CMB observations at / < lO"' will only constrain the combination S^kj. To then convert the 
constraint to a bound on the energy density in magnetic fields requires knowledge of the exponent n' and the inertial 
and dissipation scales. In other words, the CMB signature for I < 10'* probes the long wavelength tail of the magnetic 
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spectrum and not the modes where the bulk of the energy density resides. This suggests that it may be favorable to 
investigate the CMB at yet higher indeed, our results do show stronger signatures with growing I. 

Big bang nucleosynthesis (BBN) is sensitive to the total energy density in the magnetic field since this is what enters 
the expansion rate of the universe. The constraint from BBN [29-32] is best expressed in terms of the magnetic field 
energy density relative to the photon energy density 

nBy^ — < 10^1 (9) 
Pi 

The relative density, Qb-j, is independent of cosmological epoch since both magnetic field and photon energy density 
scale as where a{t) is the cosmological scale factor. Note that is the photon density which is different from the 
radiation density at BBN epoch since neutrinos also contribute to radiation. 

Finally, all the scaling due to the expansion of the universe can be pulled out by converting to rescaled quantities 

Ba^ — > B , p^a"* — > , fca — > fc , (10) 

and, in what follows, we will use these comoving quantities unless explicitly stated. The final form of the magnetic 
field power spectrum in the inertial range is 

S{k) = nB^p^^-^i^-j , k<ki (11) 

The power spectrum S{k) will enter the calculation of C^^ in combination with powers of k. So it is convenient to 
introduce the dimensionless "power spectrum" using powers of k and also the wavelength of observed radiation, Aq, 

2n 

^ I A? ( A 



167r2e j " 1 ^0 (ff) kj < k < k^i, 



^i,ik)^k^S{k){^] ={.2fj^Y-' (12) 



where 



^0 k > fcdii 

Qt) 

At the present epoch p^{to) = 4.64 x IQ-^^gm/cm^ 2 x lQ-^^{eYf [33], and so 

where we denote the observed CMB frequency by vq. Note that Aq is independent of fc/. 

In principle, fcdiss is not an independent parameter. One can estimate its value for a given amplitude and shape 
of the magnetic fields spectrum. According to [34, 35], fcdiss is determined by damping into Alfven waves and can be 
related to Beff as 

fcdiss -1.4 ^1/2 AO-^Gauss \ ^^^^ 



Converting this to ^B-f, we obtain 



/ 10^2/, 

fcdiss « 0.43 J— Mpc-i , (16) 

V "S7 



where it was useful to know that 1 Gauss = 6.9 x 10 GcV^ when converting between natural units and CCS units. 
We also note that Eq. (4) implies: 

Boff = 3.25 X 10"^ Gauss . (17) 

The relation (15) is based on the analysis in Ref. [34] where small perturbations on top of a homogeneous magnetic 
field were treated. To extend this analysis to a stochastic magnetic field with little power on long wavelengths, Ref. [35] 
introduced a smoothing procedure and split the spectrum into a "homogeneous" part and a "perturbations" part. 
It is not clear to us if this procedure is valid for an arbitrary spectrum, S'(fc), but we will still use Eq. (15) as an 
approximate expression for the dissipation scale. 
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III. FARADAY ROTATION CORRELATORS 



The CMB is linearly polarized and an intervening magnetic field will rotate the polarization vector at a rate given 
by: 

da = X^- ^aUeB-dl, (18) 

znm^ 

where A is the wavelength of light, a is the scale factor normalized so that Otoday = 1, ne is the number density of free 
electrons, dl is the comoving length element along the photon trajectory from the source to the observer and we are 
using Gaussian natural units with H = c ~ 1. Using the known expression for Thomson scattering cross-section, 

and integrating along the line of sight, we obtain the Faraday rotation of the polarization angle, 

3 .2 



^ 2 A / f(x)B.dl (20) 

where 'r(x) = UeCrTa, is the differential optical depth along the line of sight, Aq is the observed wavelength of the 
radiation and B = Ba^ is the "comoving" magnetic field. The limits of the integral are from the initial position of 
the photon to the final position. 

FR depends on the free electron density, which becomes negligible towards the end of recombination. Therefore, the 
bulk of the rotation is produced during a relatively brief period of time when the electron density is sufficiently low 
for polarization to be produced and yet sufficiently high for the FR to occur. The average FR (in radians) between 
Thomson scatterings due to a tangled magnetic field was calculated in Ref. [36] and is given by 

^=i5»o.oar^)(^r, ,21, 



87r2e vl ■ VlO^^G / V ^^o 

where is the current amplitude of the field and I'o is the radiation frequency observed today. 

In this section we will calculate the two-point correlation functions of the FR angle, a, and this will be related to the 
two-point correlation function of the magnetic field as given in Eq. (1). The FR correlator will enter the calculation 
of Cj^^ in the thin LSS approximation of Sec. IV, in which one assumes that all of the polarization was generated 
at once independently from the FR. In the general case, which we present in Sec. V, the generation of the CMB 
polarization and its FR are entangled. 

FR is sensitive only to the line of sight component of the magnetic field, whereas magnetic helicity, described by 
the helical power spectrum A{k) in Eq. (1), depends on all 3 components of the magnetic field. So a correlator of FR 
cannot be sensitive to the helical properties [22, 37, 38] ^ . Following Ref. [39] without the helical term, we get 

(a(n)a(n')) = f ^^S{k) [ dj] [ dr/f{j])T{7]')e-'^"''e+'^"'''' [fi ■ n' - {k ■ f^^ (22) 



IGTT^ey J (27r 

where n, fi' are two directions on the sky. Note that, as is conventional, we have written 'r(x) in Eq. (20) as t(?7) 
in the integrand of Eq. (22) even though x = n(?7o ~ '7)- We have also ignored inhomogeneities in the free electron 
density along different directions on the sky since this will only give a higher order correction to FR. The limits of the 
1], rj' integrations are from the time that the photon last scattered to the present time. In general, the last scattering 
time will be different for different photons but, in the thin LSS approximation, the initial time will be taken to be rj^,, 
the epoch at which the "visibility function", g{r]) = re^'^, is maximum. 

Statistical isotropy implies that {a{h)a{h')) must be a function of h ■ n' . This can also be seen directly by writing 
Eq. (22) as 

(a(n)a(n')) = (^^) J ^S{k) J drj J d^^' T{r^)f{r^') [h ■ h' - duM j d2fce-^'' "''e+''' " , (23) 



^ Indirectly though, hehcity does enter the FR, signature because magnetic hehcity plays a crucial role in the evolution of magnetic fields 
and the exponent n in S{k) (see Eq. (f 1)). 



6 



and using 

/ d2fce-»k.n„g+»k.ft'V ^ inY.i'^L + l)jLikv)jLiW)PLin ■ n') . (24) 

where jhix) are Bessel functions and Pl{x) are Legendre polynomials. Hence, (a(n)a(n')) depends only on h ■ fi' as 
expected. As a consequence, the correlator of a's can be expanded into Legendre series 

{a{h)a{h')) ^ J2 ^^^^CrPLin ■ n') (25) 

L 

and the FR correlators can also be written as {ct*i^M^L'M') — C2'^Sll'6mm' , where ulm are the coefficients in the 
spherical harmonic decomposition of the FR angle, 

"(^i) = X! ^LMYLMin) . (26) 

L,M 



Now combining (23) and (24), and introducing j']^{x) = dxjL{x), we can write 
(a(n)a(n')) = ^ J k'dkS{k) J drj J dry'f (ry)f(77') 



(2L + 1 

L 

Applying the recursion relation 



X E [(" • "')^^(^ • f^')jLikv)jL{kv')-fLikv)j'L{kv')PLin ■ h')] . (27) 

L 



[L + 1)Pl+i{x) = (21 + l)xPLix) - LPl-i{x) (28) 
to the (n • h')PL{n ■ n') term in (27) results in 

(a(n)a(n')) = (^Sfe) ^ J k'dkS{k) J drj J dr^'r{fj)r{r^')Y,^-^^^PL{h ■ h') 

^^jL-iikv)jL-i{krf) + ^^jL+i{kv)jL+iikv') - jUMiUW) • (29) 

In analogy with the way Cl's are evaluated for CMB, it is convenient to introduce "transfer" functions TL{k) and 
T^^\k) defined as 

TL{k) = / dr] T{r))jLik{ria-T])) 

(30) 



T^'\k) EE r dry f(r;)jafc(% - r?)) 



where, as defined above, 77* is the epoch at which the visibility function is maximum and 770 is the present epoch. 
Using these transfer functions in (29) and comparing to (25) allows us to write 



2L + 1' ' " 2L + 1 



(31) 



This gives 



The function T^^' can be expressed in terms of T functions by using the relation 

Ji(^) = [LjL-i{x) ~{L + l)jL+i{x)] . (32) 

ri'^ - [LTl-1 -{L + 1)Tl+i] . (33) 

The functions Thik) are independent of the magnetic field, and can be easily evaluated numerically using the ionization 
history from CMBFAST. 

Eq. (31) is our final result for the FR correlation function. It will be useful in Sec. IV where we find Cj^^ in the 
thin LSS approximation. However, the result is not useful in the general case of a thick LSS because then Cj^^ is not 
simply related to C£". 
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IV. FARADAY ROTATION OF CMB IN THIN LSS APPROXIMATION 



In the limit of instant last scattering one assumes that all of the polarization was generated at the peak of the 
visibility function. Since we are interested specifically in the FR effects, we will neglect primordial tensor modes and 
any actively sourced vector and tensor modes (including those sourced by magnetic fields) so that only E mode is 
produced at the instant of last scattering. At subsequent times, because of the residual presence of charged particles, 
some of this E mode will be Faraday rotated into B mode. To estimate this effect, we can start with Eq. (6) of [40] 
(same as Eq. (20) of [41]) which gives the B-mode coefficients 

Blm = 2(-l)'"^ ^ a-LA/i^/amaCi^fama^its ; (34) 

LM 121712 

where Ci^famj ^-nd H^^ are defined in terms of Wigner 3 — j symbols as [41] 

.LM ^ ( ..„,/ (2? + l)(2£ + l)(2/2 + l) 7 I L h\ 

Qml27n2 - ( '^f y \-m M 1712 J ^ ' 

^^.-(20-2)' (36) 

and the summation is restricted to be only over even L -\-l2 + l. From the above, we can derive the expression relating 
Cf ^ to Cf ^: 



{Bvm'Blm) - 4^ ^ ^ ^ ^tml2rn2Bu2^fi„%^rn'^H[/i,^{alME!^rn2^L'M'Et 



LM L'M' l2m2 l'2i'^2 



L ^ h 

which assumes a statistically isotropic stochastic magnetic fields with a FR angular spectrum C£" given by Eq. (31). 
Therefore, the B-mode angular spectrum in the thin last scattering approximation is 

CP^ = - ^(2L + l)Cr E(2^2 + l)Cfi'^{Hh2f ■ (38) 

It is instructive to put the expressions for the CMB observables in a form that separates the well-established physics 
of FR of CMB polarization from the particular form of the magnetic field spectrum. For example, substituting (31) 
into (38), we can re- write the latter as 

iBS _ 2 f dk 2 



cr = - J -^Kiimm. m 

where Wi{k) are "window functions" defined as 

W,{k) = 45: (2^i±lK2^(^. fL^^^_^^^^^2 ^ X + 1 j^^^^^^^j, _ [Tt'\k)f] . (40) 



47r y ill' H \2L + V ^ " 2L + 1' 

hL \ < < / 

They describe the amount of power a given wavelength k of the magnetic field spectrum contributes to a given angular 
scale I of the B-mode polarization spectrum. We note that Eq. (39) relating to the magnetic spectrum is formally 
independent of the thin LSS approximation - the approximation is used in the calculation of the window functions. 
Transfer functions 71 (fc) can be found numerically using the differential optical depth calculated in CMBFAST [42]. 
Having evaluated the window functions once, one can store them and use (39) to calculate for different choices 
of the magnetic spectrum. 



V. EXACT (THICK LSS) TREATMENT OF FARADAY ROTATION 



The thin LSS approximation decouples the process of generation of the CMB polarization by Thomson scattering 
from its subsequent FR by magnetic fields. It effectively assumes that the background E-mode polarization on all 
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angular scales was created at a single instant in time corresponding to the peak of the visibility function. While this 
may suffice for order of magnitude estimates, polarization in different parts of the sky was created at different times 
and any choice of a single time is essentially arbitrary. Furthermore, the amount of the FR strongly depends on the 
choice of the initial instant, since the amplitude of the rotation is directly proportional on the rapidly decreasing free 
electron density. In this Section we derive the exact FR window functions, denoted by Wi{k), by solving the radiative 
transport equations for the generation and propagation of CMB polarization in the presence of FR by stochastic 
magnetic fields. The form of Eq. (39) relating to the magnetic spectrum will be the same. 

The evolution of CMB Stokes parameters is described by Chandrasekhar's radiative transport equations [43]. In 
the absence of FR, equations for the q Fourier mode of linear polarization parameters Q and U are [44] 



P± + iqnP± = -tP± + S± , 
where P± (q, n,r]) — Q ± ill, h is the direction of the line of sight, fi = q ■ h, and 

2 



(41) 



(42) 



a=-2 



In the above, ±2^im("-) are spin-2 spherical harmonics, and p('''(q, r/) — [62 ~ \/6E2 ]/10, where 82 and E2 are 
the quadrupole moments of the CMB temperature and E-mode polarization for scalar (a = 0), vector (a = ±1) and 
tensor (a=±2) modes. Assuming that polarization generated by vector and tensor sources is negligible, we have 



FR rotates Q into U, and U into Q, leading to a new term on the right hand of (41) [45] 

P± + iq^J.P± = -rP± T 2iujBP± + S± . 



(43) 



(44) 



where ^^(n, 77) = /fB(r) • n, r = (rjo — r])n, f = 3Ag/(27re), and Xq and B are the comoving wavelength and magnetic 
field strength. The ordinary differential equation (44) has the inhomogeneous solution 



(45) 



with 



= -fe-^V6P(°)(q,7y)^(-z)V47r(2; + l)[e|°)(g(r;o - v))±ipr{q{m - v))] ±2>io(n) , 



(46) 



and r = drj' f. In the above, we used the identity (Eq. (16) of [44]) 



b2l^2o(n)e^'--"'- = Y.^~i)'./^{2r+l)[ef\qr)±ipf'\qr)] ^Mh) , 



(47) 



with (in what follows we will not need /3;*'°'') 



.(0) 



(x) 



^ 3 {l + 2)\ji{x) 
~V 8(^-2)! 



For small ws, we can write (45) as 

rvo r rvo 

P±{l,n)^ / dri s±{q_,h,'i]) 1 T 2* / coBdi]' 

Jo I Jv . 



(48) 



(49) 



Next, we can use the total angular momentum formalism of [44] to derive an expression for C^^ in terms of the 
magnetic spectrum S{k). From Eq. (55) of [44] we have 



(50) 



m— — / 
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where q = z. Note that in Eq. (55) of [44] the sum over m runs only from —2 to 2 because these are the only modes 
that can be sourced by scalar, vector and tensor fluctuations in the metric. However, the FR effect on the propagation 
of photons is not via perturbations of the metric tensor. Hence, to stay general, we keep the sum to be over all m 
modes. Inverting (50) and using (49) we obtain 



X I]H)'V47r(2;i + l)[e||'^(+2ll,o[+2i^m]*+-2>^/io[-2ll™]*) / u^sdv' , (51) 

where we have assumed that FR is the only source of B-mode. The angular spectrum C^^ can be written in terms 
of B;^'"^(q) as [44] 

{21 + ifCr (S^>(q)i?(")(q)) . (52) 



Introducing g{if) = fexp(— r) and Xj™ —+2 YiiQ[+2Yim]* +-2 i1io[-2^im]*, and substituting (51) into (52) we obtain 

^ey^ di^g{rj)J^ d,y'.9(ry')(^<°)*(q,^)^<°)(q,^')) 

/ 

m— — l li I2 

dnj dn'Xli^;in)Xg{n'){J drj' J ^ dv^usiv" ,n)cuB{v"' ,n')) . (53) 

Eq. (27) for the equal time two-point correlation of rotation measure is easily generalized to the unequal time FR 
correlation case above. Namely, we have 

{p dij' p' d^"ujB{il".n)uB{^"\n')) = l J ^Al,{k)Y,(^^^p,{n-n')UL{k,v,v') , (54) 



where 

VlL{k,r,,r,') = ^^rL-i(fc,r;)rL-i(fc,,70 + ^^71+i(fc,,7)rL+i(fc,r;') -ri'^(fc,r;)ri'^(fc,?7') (55) 

and the transfer functions are the same as before, except for the range of the time integration: 

rvo 

TL{k,r]) = / d7fTiT]')jLik[r]o - r]']) 

Jrj 

ri'{k,ri) = I dri'T{r]')j'L{k[rio-ri']) . (56) 

Jr] 

The in (55) can be written as a sum of terms with separated 77 and 77' dependencies: 

3 

UL{k,v,il')^Y.''L\k,vWB\k,v') , (57) 

where 



(1) / ^ -r (2) /-^ + 1-T- (3) -r(l) tr:Q\ 

We can also relate (P'^°-'*(q, 77)P'^°-'(q, ?/)) to the primordial curvature power spectrum A^(g) via 

(q, ry)p(O) (q, r/)) = q-^A' iq)P^"> (q, r^)P^°^ (q, ry') . (59) 



Putting it all in (53) and integrating over the angular dependence of q, we obtain 



where we have defined 
and 

d\%{q,k) = / din g{v)V6P'^"\q,v)4?i<livo-v))u^L\k,v)- 
Jo 

Using the expression for the integral of a product of three spin-weighted spherical harmonics 



df,Y (n)Y (f,)Y ./(2/i + l)(2i + l)(2^2 + l) ^ /i L h\fh L .2 



47r V Jni M TO2 J \ ^o, —6 — c 

and the orthogonality property of Wigner 3-j symbols [46] , we have 

^ Z,^^f;,Z,^gL - 4^M. ^ ^^^^^ ^ if / + ;i+L = even, and otherwise 

mM 



Substituting this into (60) we can write 



Cr = l /f Ai,(fc)M^,(fc) 



with 



where I + h + L — even, and 



hL 



3 

SwL{k) = ^ f ^A^q)J2\dl^liq,k)f^ 
TT J q f— f 



The thick LSS window function in Eq. (66) is very similar to the thin LSS window function given by Eq. 
Noting that 



C^^^ = U'^A\q) 

J q Uo 



no 

di^g{i^)VQP^''\q,il)ef^{q{il^-i^)) 



we can see that 



C—1 



of the thick case becomes 



of the thin case if the function d\^^j^(q^k) defined in Eq. (62) is "factorizable" , that is if 

4,1 -> nf{k,^.) / dT? ff(ry)V6F(")(g,7?)e[°'(q(7?o - v)) ■ 
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FIG. 2: Window functions Wi{k) for different values of k plotted vs I as evaluated using the full LSS treatment (solid lines) 
and using the thin LSS approximation(dotted lines) . These window functions prescribe the way in which a given Fourier mode 
k of the stochastic magnetic field contributes to the multipole I of C^^ . For example, C^^ at I ~ 8000 is sensitive to A|,/ at 
k — 0.5, 1 Mpc~^ but not to smaller k. The left and right panels show the same functions plotted on logarithmic and linear 
axis respectively. 



Thus the thick LSS case reduces to the thin LSS case if we disregard the convolution in Eq. (62). 

The convolution in Eq. (62) will in general not be factorizable, i.e. the different q modes of the source E-mode 
polarization are not all created at a single time before FR took place. Instead, functions d['^\^{q,k) determine the 
relative amount by which a given g-mode, projecting into multipole li of , is distorted by the A:- mode of the 
magnetic field projecting onto multipole L of the FR distortion spectrum. 

To evaluate Wi{k) numerically, we modified CMBFAST to calculate sources Si-^L{k) on a grid in L and k. In order 
to accurately account for magnetic fields on scales up to a given fcmax in Mpc~^, one needs to evaluate the source up 
to Lmax ^ (10*Mpc)A:max- We included the modes up to 6000, and confirmed that it is more than sufficient for all 
k and L because of the exponential suppression of the source E-modes by the Silk damping. Time required for the 
evaluation of sources needed for Wi{k) was ^ 4000 CPU-hours, where the sources were sampled for 450fc bins. Once 
the sources are calculated and stored, the sums over li and L in (66) are quick to perform. 

Even though computing the exact window functions Wi{k) takes a non-trivial amount of CPU time, it only needs to 
be done once for a given cosmological model. We have made our window functions, evaluated using the ACDM model 
with WMAP7 best fit parameters [47], publicly available at http : //www. sfu. ca/^levon/faraday .html, along with 
a Fortran code that calculates C^^ for a given A^,j(fc). 



VI. RESULTS 



In Fig. 2 we show window functions for several values of fc as a function of multipole I for the thin (dotted lines) 
and thick (solid lines) LSS treatments discussed in previous sections. These results are independent of the spectral 
features of the magnetic field. Instead, the role of the window is to specify the extent to which a given Fourier mode 
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FIG. 3: The CMB B-mode spectrum from Faraday rotation evaluated in the case of a nearly scale-invariant magnetic spectrum 
with 2n = 2n' = 0.1, kj = 1 Mpc~^, and /cmb = 61 GHz, using the thin (blue dot) and exact (solid red) treatment of LSS. The 
value of the magnetic field energy fraction Qbj = 5 x 10"* corresponds to BcH ~ 0.73 x 10~^ Gauss. The black short-dash line 
is the input E-mode spectrum, the black dash-dot line is the contribution from inflationary gravitational waves with r = 0.1, 
while the black long-dash line is the expected contribution from gravitational lensing by large scale structure. 



k of the stochastic magnetic field contributes to the multipole / of Cf . One can see that each window has a peak at 
/ approximately given by I = lO^fc, where 10'* is roughly the distance to LSS in Mpc. One can also see the oscillations 
which come from the acoustic oscillations in the E-mode spectrum. 

Comparing the exact (thick LSS) windows with the ones obtained in the thin LSS approximation, we note that 
they have comparable shapes and amplitudes near their peaks, but differ significantly at I away from the peaks. This 
difference comes because of the assumption made in the thin LSS approximation that all of the E-mode was produced 
at the same time, so that the FR occurs at the same rate on all scales. In reality, E-mode is rotated as it is being 
produced and, since the rate of FR depends on the rapidly decreasing free electron density, E-mode scales produced 
at different times are rotated at different rates. 

Let us now focus on the B-mode spectra for specific choices of the magnetic spectrum. In Fig. 3 we consider a nearly 
scale invariant magnetic spectrum with 2n = 0.1. The black dash line shows the input E-mode, while the blue dot and 
the red solid lines show the B-mode spectra obtained using the thin and exact treatment of LSS. One can see that the 
exactly calculated spectrum favors the power near the peak at the cost of the power around it. Note that in this case 
the shape of the B-mode is essentially a copy of the E-mode spectrum, except for the lack of exponential damping 
on small scales. While the EE correlations are suppressed by the Silk damping, there is no exponential suppression of 
the FR generated small scale B-mode spectrum because the magnetic field is correlated on small scales. For 2n = 0.1, 
both thin LSS and exact spectra have the asymptotic form of at high I. More generally, the asymptotic exponent 
is 2n — 1, which can be relatively large for stochastic fields e.g. for n = 5/2, PCP^ oc at large I. 

The Silk damping dissipates the perturbations in the CMB temperature and E-mode on scales smaller than 9 Mpc 
or so. In our formalism, it comes through the exponential suppression of the source function P'^°\k,T]) in Eq. (43). 
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FIG. 4; The CMB B-mode spectrum from Faraday rotation in the 2n = 5, 2n' = 3 case motivated by causally generated fields, 
with ki — 1 Mpc~^ and /cmb ~ 61 GHz. The thin blue dot line shows the thin LSS calculation, while the exact (thick) LSS 
calculation is shown with solid red. The value of the magnetic field energy fraction, Q,bj = 10""^, corresponds to BeS ~ 10~^ 
Gauss. The other lines are the same as in Fig. 3. 



However, the dissipation scale of (weak) magnetic fields at decoupling time is typically much smaller [13, 34] and 
for k < IMpc"^ the magnetic field can be treated as a stiff source, i.e. we can safely assume that its evolution is 
independent of the perturbations in the photon-baryon fluid. Then, the only damping of the FR induced B-mode 
power is due to averaging over many random rotations along the line of sight. The functional form of this suppression 
is a power law and the exponent can be estimated by observing that a random superposition of N perturbations along 
the line of sight leads to a statistical reduction in the amplitude of the observed anisotropy by a factor 1 / \/N. For 
wavenumber k we have N oc Arjk, where A77 is the period of time, comparable to the thickness of LSS, during which 
FR is efficient. Thus, the power spectrum on small scales, which is the square of the FR amplitude, is suppressed by 
1/fc, which translates into the l/l suppression of the angular spectrum. 

In Fig. 4 we show the plot of the B-mode spectrum for the magnetic spectrum index 2n = 5 expected for causally 
generated magnetic fields [13]. In the figure we used kj = 1 Mpc~^. For comparison, we show the B-mode from weak 
lensing and from gravity waves, as well as the E-mode spectrum. Note that at small angular scales (high /), the FR 
produced B-mode can dominate the signal as it keeps growing as 

At present, B-mode have not been detected; there are only upper bounds. Still, even these weak bounds can produce 
constraints on the magnetic field fraction ^bj, and furthermore, the bounds will improve rapidly as CMB observations 
are made on smaller angular scales. We derive current constraints on the magnetic field energy density, flB^y, from the 
WMAP 7-year polarization data by comparing magnetic field induced theoretical CMB B-mode power spectrum Cf^ 
as given by Eq. (65), with the WMAP observed B-mode power spectrum using the statistics. We consider three 
WMAP frequency bands Q, V, and W corresponding to frequencies 41 GHz, 61 GHz, and 94 GHz respectively. The 
lower frequency bands, K and Ka, are foreground dominated and we do not include those in our analysis. We combine 
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FIG. 5: Constraints from WMAP-7 year data on the magnetic field density Qs-y, or analogously magnetic field efTective 
amplitude B^g as defined by Eqs. (4) and (17), as a function of the inertial scale wavevector kj. We consider two choices of 
magnetic spectrum, the case of a nearly scale-invariant spectrum 2n = 2n' = 0.1 (dash-dot red line) and causal magnetic fields 
with 2n = 5, 2n' = 3 (dashed red line). The BBN constraints on the magnetic field density f2_g.^ < 0.1 is shown by the dotted 
blue line. At large ki the bound is set by the theoretical relation between ^Is-y and fcdiss > ki given by Eq. (1.5) (black solid 
line). 



different frequency channels directly when evaluating the and, when evaluating the likelihood, we restrict the 
maximum value of fi^^ to the one given by Eq (16) for kj = kdiss- We used WMAP data for ^ > 32, above which the 
errors for individual ^'s can be treated as uncorrelated. The maximum multipole considered in the analysis is ^ = 700. 
For our analysis we considered two choices of theoretical magnetic spectral indices (n, n'), one corresponding to nearly 
scale-invariant spectrum 2n = 2n' = 0.1, and the other corresponding to causal magnetic fields 2n = 5, 2n' — 3. Our 
analysis assumes that the magnetic field is the only source of the B-mode signal, and ignores the possibility of other 
sources of B-modes such as inflationary gravitational waves [48], weak gravitational lensing of the CMB [49], and 
several other distortions of primary CMB along the line of sight (see Refs. [50, 51] for examples of such distortions). 

To obtain bounds on the magnetic field, we note that the CMB FR signature constrains Ag defined in Eq. (14), 
which is proportional to Qs-y/n. Now k is defined in Eq. (7) and depends on the dissipation scale, k^iss, which can 
also be related to by using Eq. (16). Therefore for given values of n and rt', k is a function of the inertial scale, 
kj, and the energy fraction in the magnetic field, flsf In Fig. 5 we present the 95% confidence level (C.L.). bounds 
on for two choices of (n, n'): the nearly scale-invariant case, 2n = 2n' = 0.1, and the causal case, 2n = 5, 2n' ~ 3. 
The bound in each case is presented as a function of kj, which we treat as an unknown parameter in the model. 

The dependence of the constraint on ^Is-y on kj in each case can be readily understood. In the 2n = 5, 2n' = 3 
case, is independent of kj at very small kj because for positive n the integral in Eq. (65) is dominated by the 

contribution from k ^ kj. The dependence of on kj becomes stronger when kj gets closer to the range of scales 
with significant E-mode power (fc ^ 0.1 Mpc~^), which is encoded in the shapes of the window functions. However, 
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this happens to be just under the maximum wavenumber constrained by WMAP's measurement of BB. Thus, WMAP 
is unable to probe the dependence on fc/, and the constraint hne in Fig. 5 stays almost horizontal even at kj ~ 0.1. 

— 1/2 

The condition kj < fc^iss oc imposes its own upper bound on i^B-y at large kj. Namely, at sufficiently large fc/, 

the constraint curve becomes independent of the CMB constraints, and is a consequence of the constraint arising from 
dissipation. The curve kj = k^iss is also shown in Fig. 5. 

In the scale-invariant limit, which is the 2n = 2n' = 0.1 case in Fig. 5, the CMB constraint is independent of kj 
because all k dependence effectively disappears. At large kj, the bound is eventually dominated by the fc/ = fcdiss 
curve, which is independent of the CMB data. 

We note that a CMB experiment which can measure B-modes aX £ ^ 1000 would be sensitive to changes in kj in 
the case of a causal spectrum. We will present forecasted bounds from future CMB data in an upcoming paper [52]. 

VII. CONCLUSIONS 

Primordial stochastic magnetic fields may be produced in the early universe during baryogenesis. Characterizing 
such primordial magnetic field is extremely valuable for probing early universe physics. In this paper we have calculated 
the effect of primordial magnetic fields on the CMB polarization. This study is especially timely in view of upcoming 
and next generation of CMB observations that are focused on measuring the polarization of the CMB. 

FR of existing E-modes of the CMB can generate parity odd B-modes whose amplitude will depend on the obser- 
vation frequency bands. We generalized the signatures previously obtained assuming "thin" last scattering surface by 
solving the full CMB radiative transport equation with FR taken into account. Comparing our full treatment with the 
thin last scattering approximation, wc find the respective B-mode spectra have similar shapes but their magnitudes 
at a given multipole £ may differ by as much as a factor of a few depending on the primordial magnetic field spectrum 
and the angular scale of observation. The well-established physics of FR of CMB can be encoded in window functions, 
Wi{k), which determine the relative amount that a given Fourier mode k of the magnetic field contributes to the 
multipole I of C^^ . The window functions are independent of the details of the magnetic field spectrum and only 
need to be computed once for a given cosmological model. We have evaluated them in the best fit ACDM model and 
made them available at http://www.sfu.ca/^levon/faraday.html, along with a Fortran code that calculates C^^ 
for a given magnetic spectrum. 

FR can take a smooth E-mode field and distort it to create B-modes on arbitrarily small scales. For instance, on 
scales smaller than the Silk damping scale, the E-mode map is essentially homogeneous. However, FR creates E and 
B-mode inhomogeneities on sub-Silk scales by rotating different parts of the homogeneous E-mode patch in different 
random ways. As explained in the previous section, the power of these inhomogeneities is suppressed as l/l due to 
random superposition of multiple small scale rotations along the same line of sight. 

FR of the CMB due to magnetic fields has a distinct frequency dependence z^^) which will allow it to be 
distinguished from other sources of CMB B-modes such as lensing, inflationary gravitational waves, topological defects 
[53, 54] and FR due to pseudo scalar fields [55-58]. The growth PCj^^ oc at large I is also characteristic of a 

magnetic field and will help discriminate primordial FR from other foreground contamination. We have summarized 
current constraints on magnetic enegy 0,3-^ from WMAP7 based on the FR induced B-mode spectrum in Fig. 5. For 
comparison we have also shown the constraints on magnetic energy obtained by BBN. For a scale invariant magnetic 
field, the constraints from WMAP 7-year data are 2-orders of magnitude better than the BBN constraints. For a 
causal magnetic field with (2n ~ 5, 2n' = 3), BBN constraints are weaker only by a factor of few. However, note that 
the constraints from BBN are not expected to improve much in future, but future observations of CMB polarization 
at smaller angular scales will significantly improve the constraints on magnetic field. We also refer the reader to [59] 
for the most recent CMB bounds on magnetic fields based signatures other than FR. 

Although in this paper we have focused on B-mode power spectrum of CMB generated due to FR of the CMB 
polarization, recently it has been shown that the stochastic FR of polarization of the CMB couples different CMB 
angular modes, thus generating non-Gaussianity which can be seen in the CMB trispectrum [40, 41, 60]. In a follow-up 
paper [52] we discuss the detectability of primordial magnetic field using such non-Gaussian features in the CMB. 
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